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We follow and modify the Feshbach-Villars formalism by separating the Klein-Gordon equation 
into two coupled time-dependent Schrodinger equations for particle and antiparticle wave function 
components with positive probability densities. We find that the equation of motion for the proba- 
bility densities is in the form of relativistic hydrodynamics where various forces have their classical 
counterparts, with the additional element of the quantum stress tensor that depends on the deriva- 
tives of the amplitude of the wave function. We derive the equation of motion for the Wigner 
' function and we find that its approximate classical weak-field limit coincides with the equation of 

£\j , motion for the distribution function in the collisionless kinetic theory. 

o 

Q 

' I. INTRODUCTION 

(N ' 

In relativistic wave mechanics, it is well known that the naive probability density p = 2Im(?/>* dtip) , as constructed 
from the wave function ip of the Klein-Gordon equation, is not necessarily a positive quantity. The presence of a 
Oh' negative probability density may appear to preclude its description in relativistic hydrodynamics. The origin of such 
. pathology arises because the wave equation is second order in the time-derivative, and the corresponding probability 
' density involves a time-derivative of the wave function. Dirac's resolution of the problem led to the introduction of 
another wave equation, the Dirac equation, that is first order in the time-derivative [H, HI • A parallel resolution of 
the pathology was provided by Pauli and Weisskopf in quantum field theory with the introduction of particles and 
antiparticles, and the wave field ip is interpreted not as the probability amplitude but as comprising of operators that 
create and destroy particles in various modes Q. 

Wave mechanical interpretation of the wave field tp was however reintroduced by Feshbach and Villars Q by noting 
that the Klein-Gordon equation actually contains both particle and antiparticle degrees of freedom. The particlc- 
' antiparticle separation can be achieved by writing the Klein-Gordon and the Dirac equations as a set of coupled time- 
dependent Schrodinger equations for the particle and antiparticle wave function components with positive probability 
densities 

In many practical problems, as for example in the evolution of dense matter with relativistic constituents in 
relativistic hydrodynamics [oT-floj or in the quarkonium two-body problem (ill Il2j , particles and antiparticles can be 
approximately treated as two distinct types of interacting particles, and they possess positive probability densities. 
For these problems, we are motivated to follow the Feshbach-Villars formalism where the particle and antiparticle 
J^j ' probability densities can be positive definite. 

We wish however to modify the formulation of Feshbach and Villars Q • In their formulation, the relativistic prop- 
el ■ erties of the momentum variables in the coupled Schrodinger equations are not readily apparent and the Schrodinger 
equation for the antiparticle wave function contains the operator — (p — eA) 2 /2m that differs from the standard kinetic 
energy operator by a sign. It will be desirable to reformulate the problem so that the relativistic properties of the 
momentum variables become more apparent and the Schrodinger equations for the particle and antiparticle wave func- 
tions contain kinetic energy operators with the same sign. As a result, the connection to relativistic hydrodynamics 
can be more readily worked out. 

Using the new set of coupled Schrodinger equations, we wish to search for a hydrodynamical description for the 
evolution of relativistic probability densities. Hydrodynamics and quantum mechanics have many elements in common, 
since the density fields and the velocity fields are important dynamical variables in both descriptions. It is therefore 
well-known that the Schrodinger equation can be cast into a hydrodynamical form for the evolution of the probability 
density [HGJ]. Such a correspondence has been utilized [l5[ to form the foundation for an "a posteriori' theoretical 
support for the validity of treating a nucleus as a liquid drop, as in Bohr and Wheeler (T(| . and treating the fission 
of a nucleus in liquid-drop hydrodynamics, as in Hill and Wheeler [l7j . There are however important differences 
associated with quantum effects that are absent in classical hydrodynamics [3, HH, [H, [3 ■ The deviation of the 
classical hydrodynamical description from the quantum treatment is embodied in the presence of the quantum stress 
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tensor pf^ in quantum fluids [l4l . [l5j . Quantum shell effects manifest themselves as nuclear shell effects superimposed 
on a smooth hydrodynamical liquid-drop background and they lead to the intrinsic deformation in many nuclei fl9j . 

In the related area of hadron and nuclear collisions, a relativistic hydrodynamical description of the collision process 
is a reasonable concept, as pioneered by the work of Landau [f| and supported by recent experimental findings |20l - [22l |. 
Relativistic hydrodynamics has been applied to study the evolution of a quark-gluon matter in the work of Bjorkcn 
Q, Baym et al. 0], Ollitraut @, and many others @. They led to successful investigations on the dynamics of 
matter with relativistic constituents in extreme conditions, as occurs in relativistic heavy-ion collisions [2(], Hil- Wave 
mechanical description of quantum systems in terms of probability densities and probability currents contains the 
proper theoretical ingredients appropriate for continuum hydrodynamics. It will therefore be of interest to generalize 
previous formulation of hydrodynamics in the Schrodinger equation in Ref. [l4| to Klein-Gordon and Dirac equations, 
in order to investigate how the wave mechanics of Klein-Gordon and Dirac particles may be used to provide the 
foundation for relativistic hydrodynamics of dense and compressed systems with relativistic constituents. 

As the foundation of hydrodynamics is usually presented within the framework of the kinetic theory [23 - |27j , it is 
instructive to investigate how the present quantum probability density approach and the kinetic theory approach are 
connected. We shall study how the equation of motion for the Wigner function, derived from the time-dependent 
Schrodinger equations representing the Klein-Gordon equation, can be related to the equation of motion for the phase 
space distribution function in the kinetic theory, in the classical weak-field and collisionless limit. 

This paper is organized as follows. In Section II, we introduce the particle and antiparticle wave functions to 
represent the Klein-Gordon wave function and its time derivative. The Klein-Gordon equation is then separated into 
two coupled time-dependent Schrodinger equations. The particle and antiparticle probability densities obey equations 
of continuity containing additional terms, but the total net particle number remains a conserved quantity. In Section 

III, we examine the Euler equation for the motion of the probability fluid of a particle or an antiparticle. To provide 
insight into the dynamics, we specialize to the simplified case when the particle-antiparticle pair production can be 
suppressed. The motion of the probability fluids of particle and antiparticle obeys relativistic fluid dynamics equations, 
and the quantum stress tensor provides part of the source of the equation of state of relativistic matter. In Section 

IV, we examine the generalization to relativistic hydrodynamics for a simple many-body system in the mean-field 
description. In Section V, we derive the equation of motion for the Wigner function and compare it with the equation 
of motion for the distribution function in the kinetic theory. In Section VI, we present our summary and discussions. 



II. SEPARATION OF KLEIN-GORDON EQUATION INTO PARTICLE AND ANTIPARTICLE 

COMPONENTS 

In the present investigation, we are interested in the temporal evolution of a single-particle boson state v with an 
initial wave function tp{ r ^) at t = 0, in external fields that consist of a scalar field S and a vector gauge field (A , A) 
with a coupling constant e. The single-particle state is characterized by a (net) particle number n v that is a conserved 
quantity (see Eqs. (|1 5|) and (|16[) below). The particle number of a state is quantized; it has the value of n v = 1 when 
v is a particle state and n v — — 1 when v is an antiparticle state. The evolution of the state v is described by the 
Klein-Gordon equation 

(ihdt - eA ) 2 ^ = [(-V - eA) 2 + (m + S) 2 ]^. (1) 
i 

To separate out the particle and antiparticle degrees of freedom, we introduce the energy function E and the auxiliary 
wave function ip^. For a given wave function ip for the state v with a charge e v = n„e, the energy function E is the 
positive root of the following quadratic equation of E, 

J dr r{(E - e v A°f + [ihd t {E - e v A )] - [(^V - e,A) 2 + (m + S) 2 }}^ = 0, (2) 

which is obtained by taking the scalar product of the wave function with the Klein-Gordon equation (TTJ. We may 
not know dtE initially at t = 0, but E and dtE can presumably be evaluated self-consistcntly and iteratively, when 
we succeed in getting the equations (Eq. ([7]) below) that allow us to advance the wave function to the next time step. 
Following Feshbach and Villars we introduce the auxiliary wave function V>4 to represent the time-derivative of 
the wave function ip, 

(itidt - eA°)4> = (E — e„A°)iP 4 . (3) 

Our definition of ip^ m the above equation differs from that of Feshbach and Villars [H, where the right-hand side 
is given as mip^. The present definition exhibits better the relativistic properties of the momentum variables and 
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facilitates the representation of the Klein-Gordon equation in rclativistic hydrodynamics. The Klein-Gordon equation 
(P) becomes 

(ihd t - eA°)i> 4 = E _\ A0 |(^V - eA) 2 + (m + S) 2 ]^ - [iW t (E - eA )}^ . (4) 
We define the particle and antiparticlc components \± of the state v as linear combinations of ip and tp4, 

x+ = + M (5) 

X- = -j=[r-1>l] or x *_ = ±^-^]. (6) 

Our definition of X- m Eq. © contains a complex conjugation that is different from the definition of X- of Fcshbach 
and Villars Q. Such a modification allows one to obtain a Schrddinger equation for x- that contains the kinetic 
energy operator (p — eA) 2 /2(E — e v A°) 1 instead of a kinetic energy operator with the opposite sign in Feshbach and 
Villars In terms of the wave functions x+ an d X- as defined in Eqs. ^ and ©, equations ^ and @ become 

{ihdt - e ± A°)x± = 2{E ^ A0) |(yV - e ± A) 2 + (m + S) 2 + [{E - e v A a f - ihd t {E - e„A°)]j X ± 

+ a^, 1 ^) {(t v - e ± A ) 2 + ( m + 5 ) 2 - ^ - e ^ ) 2 - m ^ E e ^°4 (7) 

where e± = ±e. Thus we obtain the central result that the Klein-Gordon equation can be reduced to a set of two 
time-dependent Schrodingcr equations in which the particle and the antiparticlc appear as distinct types of particles 
with opposite charges e± with the particle wave function X + an d the antiparticle wave function x~- A general solution 
of the Klein-Gordon equation contains two components that can be represented by a column vector, 

X+ 



Both x+ an d X- have positive norms, |x±| 2 i which can be interpreted as the probability densities of particles and 
antiparticles respectively, as in hydrodynamics. 

Note that our denominator in Eqs. (0 differs from the denominator of m in the formulation of Feshbach and Villars 
The relativistic properties of the momentum variable are more apparent in Eq. (|7|)- These relativistic properties 
will facilitate the subsequent representation of the Klein-Gordon equation in relativistic hydrodynamical form. 

The second term of the new equation ([7]) represents the particle-antiparticlc coupling and pair production. It is 
proportional to {{h/i)V — ej-A) 2 + (m + S) 2 — [{E — e u A°) 2 — ihd t {E — e„A )] and involves essentially the deviation 
of E 2 from p 2 + m 2 that increases with the strength of the interaction. The coupling is inversely proportional to the 
energy (or mass) of the particle. The particle-antiparticlc coupling and the rate of pair production are large, when 
the strength of the interaction relative to the rest mass of the particle is large [13, [H, ■ The particle-antiparticle 
coupling is small when the strength of the interaction relative to the energy (or mass) of the particle is small. 

For stationary states x± with e v = n v e = ±e in static external fields, Eq. (|7|) gives 

(E - e ± A ) 2 - ( Jv - e±A) 2 - (m + 5) 2 | X ± = (~(E- e ± A a ) 2 + (^V - e ± A) 2 + (m + 5) 2 | X * T - (9) 

If the right-hand side representing the pair production coupling can be neglected, the above equation is just the 
time-independent Klein-Gordon equation for a stationary state of a particle (or an antiparticle) in external fields. 
Equation ([7]) has the correct limit for stationary states in static external fields in the absence of pair productions. 

In a static case with time- independent external fields, the energy function E is a constant of motion, although 
the wave functions X ± an d the associated probability fluids in the static fields can still evolve dynamically in space- 
time. With time-dependent external fields, the energy of the single-particle state changes because energy may be 
supplied or removed by the external fields. We envisage conceptually that a two-component wave function ^(r, t) = 
(x_l_(r, i),x_(r, t)) for the state v with e„ = n v e is initially known as t = 0. The knowledge of the initial wave 
functions allows us to obtain E from Eq. ([2]) and to construct x+ an d X- a t the next time step with the help of Eq. 
([5]) and In the initial evaluation of E and in the subsequent stepwise increment of the wave function in time, we 
need the quantity d t E(t), which may be self-consistently determined by an iterative procedure. 

After the wave function ^> has been advanced to the next time step, we can evaluate E{t) using Eq. (0 that is a 
quadratic equation in E. It contains both a positive- E and a negative-i? solution. We shall limit our attention only 
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on the positive-i? solution so that we can speak of both particle and antiparticles states of positive energies. This 
stepwise increment allows the determination of the evolution of the state \&(r, t) and E(t) as a function of time. 

To see how the probability fluids behave in space and time, we write the wave functions x ± in terms of the amplitude 
and phase functions, 

X±(r,t) = ^±{r,t)e lS ^ r ' t ^ R - in ^. (10) 
We construct X± x ©~X± x ©*• After some manipulations, we find 

d t [(E - e u A°)4> 2 ± ] + V • [4>l(VS ± - e±A)] = X±, (11) 

where 

2A± = M(^V- e± A) 2 X ;-x±(^V-e ± A) 2 XT } 
+ [(m + S) 2 + (E- e v Af]{ X ±X% - X±X T ) 

+ [ihd t {E - e v A°)](x±xl + X±Xt)- (12) 

As the quantities X + and X_ are not generally a full divergence, the total number of particles and antiparticles in 
the two components are not conserved. However, the difference of the particle number and antiparticle numbers of 
the two components satisfies the equation 

d t [(E - e„A°)(<j? + - 4> 2 _)] + V • [<fi.(yS+ - e+A)] - V • [0 2 _(VS_ - e_A)] = X+ - X_. (13) 

But, X + -X^ is a complete divergence, 



X+ - X_ = - V • (x;V X *_ - X *-V X ;)/2 + V • (x+V X - - X-V X+ )/2 
Therefore, the quantity 



V-[e + A°( X * + x*-+X+X-)/i]- (14) 



E-e v A° /l2 ... [E-e u A\, 



n part icic = / dr ^— {<t>\ - r_) = / dr — (|x+| 2 - Ix-H (15) 

J m J m 

is a conserved quantity because the right-hand side of the equation (TTJJ is a complete divergence. The additional 
number of particles produced is equal to the additional number of antiparticles produced. The equal increase in 
particle and antiparticle numbers associated with X + and X- represents the occurrence of particle-antiparticle pair 
production. 

A single-particle solution with a n par ti c ie = n v = ±1 is one in which |x±| 2 3> |XtI 2 anc l can l* e normalized to be 
f E — e„A° 

/ d 3 r — [|x±| 2 — Ix=f| 2 ] — 1 f° r a particle state with n par ti c ic = ±1- (16) 

J m 

III. KLEIN-GORDON EQUATION IN HYDRO DYNAMICAL FORM 

To obtain the Euler equation, we construct X± x 0+X± x 0*- Putting all terms together, we get 



b 2 ± (-2d t S± - 2e±A°) 
1 



[20 ± V 2 </> ± -20 2 t (V5 ± -e±A) 2 ] 



2(E - e v A°) 
+[(m + S) 2 + {E - e v Ay]24> 2 ± 

+Xi(-V - e±A) 2 x; + X±(— V - e±A) 2 x T 

i —i 

+ [(m + S) 2 -{E- e„A) 2 }(x* ± X* T + X±X T ) 



+[ihd t (E-e u A°)](x*±X% - X±X T ) ■ (17) 
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In the present hydrodynamical description, we investigate the evolution of the probability densities of particles 
and antiparticlcs in a single-particle system with a definite particle number n v in external fields. The rate of the 
particlc-antiparticle pair production depends on the strength of the external interaction relative to the particle rest 
mass. In the case of a strong interaction and large pair production probabilities, a hydrodynamical description will 
be complicated as it will involve a component that describes the pair production process. 

A hydrodynamical description will be appropriate after the active pair production stage has passed and the expan- 
sion of the system is driven by a slowly varying external field or a mean field. It is this type of motion for which we 
wish to provide a hydrodynamical description. This is the case when the interaction is present but not large com- 
pared to the rest mass, so that the probability of pair production is small, and can be suppressed in the lowest-order 
approximation. This is equivalent to the case of a "simple fluid" in relativistic hydrodynamics, in which the chemical 
composition of the fluid ceases to change [30| . In such circumstances, one can speak of a single-particle system with 
a definite particle number n v , which can take on n v = 1 , — 1 values. We shall now consider this simplifying case with 
suppressed pair production by discarding the last four terms in the curly bracket on the right hand-side of Eq. (jTTJ) . 
After dividing by —2§\, the equation for the phase function S± for this simplified case is 



(d t S ± + e ± A°) 



1 



2(E-e„A°) 



[(V 2 </>±)/</)± - (V5± - e±A) 2 ] 



-(m + S) 2 - {E~e v A° 



(18) 



For this case with suppressed pair production, e v 
and multiply by 4>±(E — e±A°). We obtain 



n v e = e±. We take the gradient Vi of the above for i = 1,2, 3, 



(E-e±A )<t> 2 ± d t (V i S±-e±A i ) 



^ 3=1 



-(m + S)ViS - 0±(Vj£± - e±A*)e±F i - 

3=1 



(E - e±A°)<f>le±F oi 



e±V 1 AV| 
~2{E-e±A°) 



(V 2 0±)/0± - (V5± - e±A) 2 ] - (to + S) 2 + (E - e ± A ) 2 



(19) 



Using the equation of continuity for this simplified case without pair production, we obtain 



(_E- e± A>|(V,5 ± -e±i' 



TO + S 



E^ 

3 = 1 



%(V 3 -S± - e±A'){VjS± - e±A l ) 

TO + S 



TO + S 



3 3 

E V ^S } - + - e±A°)0 2 ± e ± F 01 - ^ </> 2 ± (V,S ± - e ± A^)e ± F^\ 

3=1 ^ 3=1 ' 



e±ViA°<j) 2 ± 



2(E - e ± A°)(m + S) 



(V 2 </)±)/^± - (VS± - e±A) 2 ] - (to + S) 2 + (E — e±A°f 



-{E-e^miV^-e^ 1 )- 



dtS 



(to + S) 2 

We can identity the fluid energy density e± as 



3 

E 

3=1 L 



ir±(VjS± - e±Ai){ViS± - e±A l 



m + S 



Vi«S 



(to + Sf 



(20) 



e± = {m + S)4>±, 



(21) 



as it corresponds to the fluid energy density for the fluid clement at rest. The fluid element is characterized by a 
relativistic 4- velocity u M . We can identify 



u± 



E 



.4° 



TO + S 

V t S± - e±A l 

TO + S 



for i = 1,2, 3, 



(22) 



G 



which obeys (u!]_) 2 — (u±) 2 = 1 in the absence of the quantum effects. We can then write an equation of motion for 
the probability densities in terms of the hydrodynamical equation 



i=i 



„(«) 



J'=l 



+ — Lj-(£ - e ± A°)</> 2 ± e ± F w - V <j>l(VjS± - e ± A j )e ± F^\ 
m + b y ' J 

.7 = 1 



e±ViA»<f>i 



2{E - e ± A°)(m + S) 



(E - e±A°)cj>l(ViS± - e±A l ) 



(V 2 0±)/^± - (V5± - e±A) 2 ] - (m + S) 2 + (E - e±A ) 2 



d t S 



3 YA2 



(m + S) 2 



E 



^(Vj-Si - e±A>)(yiS± - e±^ 
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(m + 5) 2 ' 



(23) 



where i,j = 1, 2, 3. This is the Klein-Gordon equation for the particle and antiparticle probability densities in hydro- 
dynamical form. The first two terms on the left-hand side correspond to d^T± , with the energy momentum tensor 

of the probability fluid T± — e±u±u±, for \x = 0, 1, 2, 3. The third term on the left-hand side is the quantum stress 
tensor arising from the spatial variation of the amplitude of the single-particle wave function 



14] 



„(«) 
P±ij 



— v, 

m 



>±Vj(p±. 



(24) 



In Eq. (|2"3")l there is no thermal stress tensor contribution in Eq. ([2"3"| for a single-particle. As indicated in [3], a 
thermal stress tensor pf^ will however arise from different velocity fields when there are many particles in a many-body 
system. The first two terms on the right-hand side contain forces coming from the scalar interaction, the electric field 
F° l and the magnetic field F lJ , as in classical formulations. The third term on the right-hand side is the relativistic 
correction to the time-like part of the vector interaction, and the last two terms represent relativistic corrections 
associated with the spatial and temporal variation of the scalar interaction. Thus the dynamics of the probability 
fluid obeys an equation similar to the hydrodynamical equation, with forces on fluid elements arising from what one 
expects in classical considerations. The additional element is the presence of the quantum stress tensor pjv that is 
proportional to h 2 and arises from the quantum nature of the fluid. 

While we have examined the hydrodynamical form when the particle-antiparticle pair production has been sup- 
pressed, the particle-antiparticle pair production can be included in future studies, at the expenses of increasing 
greatly the complexity of the simple picture we have obtained. 

IV. APPLICATION TO A MANY-BODY SYSTEM IN THE MEAN-FIELD DESCRIPTION 



The external mean field we have been studying can come internally from the single-particle state probability de nsity 
as in the Gross-Pitaevskii equation [3fT[3~4| or from the self-consistent mean- field in a many-body system [l4l HH [35|. 
A many-body system in the time-dependent mean-field description consists of a collection of independent particles 
moving in self-consistent mean- fields generated by all other particles 0, EE H3 • Each single-particle state is 
characterized by a state label a, particle type v, energy e„, and occupation number n a „. For simplicity, we consider 
the case in which the mean- field potential arises from a scalar two-body interaction v s (r\,T2) and a time-like vector 
interaction wo( r i, r 2)- We further neglect the last three terms on the right-hand side of Eq. (|2lt)) which represent 
higher-order relativistic corrections. The equation of motion for the energy density t av and velocity fields u x av for 
i = 1,2,3 and v = ±, in the single particle state a and particle type is, is then 



d t {e au u° av u l aLI ) + V 'je av u\ v i 



m + S 



m + S 

where, in the frame with the fluid element at rest, 



o=i 



E - e au A° ± , dA° 



S(r,i) = J d 3 r 2 n(r 2 , t)v s (r, r 2 ), 



(25) 



(26) 



7 



A°(r,t) = J d 3 r 2 jn+(r 2 ,t)e + + n_(r 2 ,i)e_ ju (r,r 2 ), 



(27) 



n v = y^Tt i>^n y , and n = n. + + n_. (28) 

a 

We consider a strongly interacting system in which the number of particles and antiparticles are equal so that 
n+(r 2 ) = n_(r 2 ) and n+(r 2 )e+ + n_(r 2 )e_ is zero. Then the contribution from the second term on the right-hand 
side of Eq. (|23|) is zero. Multiplying Eq. (|2~5j) by n av and summing over {a, v}, we get 

3 

dt(£2 n ^avu\ v u l av ) + v i(^Z n av t av u\ v v? av ) 

au j— 1 az> 

E^E^) + (E"-^)V^ = 0. (29) 

j—l av av 

We define the total energy density e by 

Y j n a ye av = e, (30) 



at' 



and the average 4- velocity u by 

= Y j n av e av u av /e. (31) 



We can introduce the thermal stress tensor pfj for {i,j} — 1,2,3 as the correlation of the deviations of the single- 
particle velocity fields from the average 

Y n av e av {u l av - u l ){u 3 au -u r )= p[f . (32) 

a v 

For the case with the suppression of pair production, we obtained the equation of hydrodynamics 

3 



ft(etiV) + £ fa («V +P$) + ^T5 V ^?} = °- 

.7=1 



(33) 



j 

where the total quantum stress tensor is 

h 2 . h 2 



4m ^ ^ rav LJ m 

C 

(v) 

and the pressure due to the interaction is 



—pt»\r,t)=n(r,t)V t S(r,t)=n(r,t)^- I d 3 v 2 n{v 2 ,t)v s {v,v 2 ). 



(35) 



(v) 

The mean- field stress tensor p\- can also be given as 



( )= f ayrun) [v) 



PW = {n v Qn ' -W^n^6 ih (36) 

where is the energy per particle arising from the mean-field interaction. As an illustrative example, we can 
consider a density-dependent two-body interaction 

v.(r, r 2 ) = [at + a 3 n((r + r 2 )/2)]<5(r - r 2 ). (37) 



8 



The contribution of the mean-field interaction to the stress tensor is then 

pg } = i[a 2 + 2a 3 n(r)]n%, (38) 

whose magnitude increases with the density and the strengths of the interaction. 

The quantum stress tensor pf^ and the thermal stress tensor pf^ can take on different values, depending on the 
occupation numbers n au of the single-particle states that determine the degree of thermal equilibrium of the system. 
The quantum stress tensor depends on the amplitudes of the wave functions while the thermal stress tensor depends 
on the phases of the wave functions and the deviation of the velocity fields from the mean velocities. The quantum 
stress tensor is less sensitive to the degree of thermalization as compared to the thermal stress tensor. In the time- 
dependent mean- field description, the motion of each particle state can be individually followed (35| . The occupation 
numbers n av of the single-particle states will remain unchanged, if there are no additional residual interaction between 
the single particles due to residual interactions. When particle residual interactions are allowed as in the extended 
time-dependent mean-field approximation (36j , the occupation numbers will change and will approach an equilibrium 
distribution as time proceeds. 

We note that the total pressure arises from many sources. We come to the observation that in situations when 
\pf) + p^\ ^ Pif for a strongly-coupled system at low and moderate temperatures, there can be situations when the 
system behaves quasi-hydrodynamically, even though the state of the system has not yet reached thermal equilibrium. 
In this case, the hydrodynamical state is maintained mainly by the quantum stress tensor and the strong mean fields. 

What we have discussed in this Section is only a theoretical framework that exhibits clearly the different sources 
of stress tensors. To study specifically the dynamics of the quark-gluon plasma, for example, it will be necessary to 
investigation the specific nature of different constituents and their interactions. Nevertheless, the general roles played 
by the different components of stress tensors can still be a useful reminder on the importance of the quantum and 
mean- field stress tensors in the strongly-coupled regime, at temperature just above the transition temperature T c . 



V. CONNECTION TO THE KINETIC THEORY 



The foundation of hydrodynamics is usually presented within the framework of kinetic theory, in which particles 
and antiparticles are described as distinct constituents and their interactions are weak [24j-[27| ■ In such a description, 
particles are considered to be approximately on-the-mass-shell, and their inter-particle collisions lead to thermalization. 
The state of the system for particles or antiparticles of type v is described by a distribution function /„(r, p;i) in 
phase space. Successive gradient expansions of small deviations from the equilibrium distribution lead to various 
approximations of the transport coefficients. The time-dependencies of the moments of various kinematic operators 
lead to hydrodynamical equations. The near-mass-shell condition restricts its application to systems with weak 
interactions that can be taken as perturbations in quantum field theory 

For a dense and strongly interacting system, such as a nucleus or a strongly-coupled quark-gluon plasma, a rea- 
sonable description of a non-equilibrium system can be formulated in a different approach, the quantum probability 
density approach considered here. In this approach, constituents of the quantum system move in the strong mean 
fields generated by all other particles. Each particle executes its single-particle motion in the time-dependent mean 
field, and the residue interaction between particles lead to "collisions" that change the single-particle state occupation 
numbers n av . These residual-interaction collisions occur in such a way that the single particle occupation numbers 
approaches a thermal distribution as a function of time [3r| . 

In the quantum probability density approach, as the particles resides in a strong field, they are off-the-mass shell 
and their energies depend on their local potential. The dynamics of each single particle state is described by a 
wave function with a positive probability density and a probability current. The equations for the total probability 
density and probability current are analogous to the relativistic hydrodynamical equations. The stress tensors in 
such a description arises from many components: the quantum stress tensor, the thermal stress tensor, and the 
mean-field interaction stress tensor. The behavior of the dynamics in such a quantum description need not be the 
same as classical hydrodynamics because the constitutive equations relating various stress components with densities 
and other attributes may respond differently to the evolving dynamics [38| . We mention earlier that in situations in 
which the quantum stress tensor and interaction stress tensor dominate over the thermal stress tensor, the degree of 
thermalization may not be important in the hydrodynamical behavior of the system. In the other extreme when the 
stress tensor arises predominantly from the thermal stress tensor, the dynamics will then depend on the degree of 
thermal equilibrium. 

Although the quantum probability density approach is particularly appropriate for cases of strongly interacting 
quantum systems, it also has a well-defined classical and weak-coupling limit that should coincide with the kinetic 
theory approach. It is therefore instructive to investigate how the two approaches are connected by studying how 
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the equation of motion for the Wigner function, derived from the time-dependent Schrodinger equations representing 
the Klein-Gordon equation, can be related to the equation of motion for the phase space distribution function in the 
kinetic theory, in the classical weak- field and collisionless limit. 

Both in the present quantum probability density approach and the kinetic theory approach, the application to 
relativistic hydrodynamics are confronted with the problem of pair production in which the interactions will lead 
to spontaneous production of particle-antiparticle pairs (28l . |29| . The standard hydrodynamics is one in which pair- 
production probability is suppressed. The equation of motion of the single-particle state of particle type v is then 
obtained from Eq. ([7J by neglecting the second term on the right-hand side. We further consider the motion to be 
sufficiently slow so that we can neglect the term d t (E — e u A°). The equation of motion for the single-particle state is 
then 



{ihd t - e u A°) X u 
which can be written as 



1 



ihdtXv 



2{E - e,A ) 



1 



2{E - e„A°) 



(-V - e„A) 2 + (m + S) 2 + [(E - e„A°) 2 } \ X u, 



( 7 V - e u A)H X u + V(r) X u, 



where 



V(r) = { + 



E 



,A° 



(39) 



(40) 



(41) 



From this time-dependent Schrodinger equation, we wish to obtain the corresponding equation of motion for the 
single-particle Wigner function. The construction of a gauge invariant Wigner function from the density matrix 
in the presence of an external gauge field (A , A) has been examined by many authors [39l - l48j , and we can follow 
similar procedures. We construct xl( r i> t)Xv( T 2, t) and introduce (ri + r2)/2 = r, and iq — r 2 = s. We define the 
three-dimensional gauge-invariant single-particle Wigner function for particle type v as 



/„(r,p;f) = [ dse^i^-'xtiv + sMxAr-s^t) 



(42) 



where the momentum p represents the kinetic momentum |43l |44| , and J 7 / is the well-known gauge factor introduced 
first by Schwinger (37j . 



pr-t-s/ z 

Ti = cxp{— ie v I A(s') • ds'}. 

Vr-s/2 

The inverse transform of Eq. (|42|) is then 
Xt(r + s/2,t) X v{r - s/2,t) 



f dn r +s/2 



(43) 



(44) 



From the Schrodinger equation (|39|) . we construct xt( r i^)*^*Xi'( r 2: t) and Xv( r 2, t)[— ihd t }xl( r i, t), subtract the 
two, and we get 



ihd t [xi{ri,t)Xv{*2,i)} 



1 



2[E - e v A°{r 2 )] 
{ 1 



-V r2 - e„A(r 2 ) 

i 



V{v 2 ) 



2\E ■ 



e v A(n) +V(r 1 ] 



xi{n,t)x»(T2,t). (45) 



The equation of motion for the Wigner function f v (r,p) can be obtained by substituting Eq. (|44|) into the above 
equation. The terms involving V(r) give the result 



V{y 2 )-V{y 1 ) 



Xt(r + s/2,t)x,(r-s/2,t) = J ^Le*»*Ti\ sin|^- ^ W(r)/ tf (r,p;t). 



(46) 



where applies only on /„ and Vjf applies only on V(r). This indicates that the quantum equation of motion for 
the Wigner function contains transcendental functions of the operators HV P ■ V r applying on the potential and the 
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Wigncr function. The expansion of the sine function will lead to a power series in h. We see here that the equation 
of motion for the distribution function in kinetic theory corresponds only the lowest order approximation in such an 
expansion. We note also that the Wigner function is in general not identical to the distribution function because it 
can take on negative values [HI, [H3|. In the classical limit, the Wigner function can be confined to be positive and be 
identified with the phase space distribution function. To make connections with the classical kinetic theory approach, 
we shall take this limit of h — > 0. The above equation is then approximated as 



V(r a )-V(ri) 



xUv + s/2,t)xu(r-s/2,t) = J -^ e -*P^[ Vr V(r)] • V p / V (r,p;t). 



(47) 



We shall take the weak- field limit so that we keep only terms first order in the external fields. Then the quantity 
V r V(r) is 



V r F(r) 



{m + S)VS (m + S) 2 e„VA° Ve^A 
(E - e v A°) + 2{E-e v A a ) 2 + 2 



+ 



(48) 



which goes to VS + Ve v A° in the non-relativistic limit of E — > m. 

To evaluate the other terms, we note that the variation of the gauge factor Ti arising from the variation of the end 
point r 2 is given by 



J Tn 



(49) 



The difference of the path integrals when one of the end points is varied from r 2 to r 2 + Sr 2 can be turned into a loop 
integral by noting that 



,[/ A(s')-ds'] = / - / 



A(s') • da' 

/T2+<5r2 

A(s') ■ da' 



A(s') • da' 
- A(r 2 ) ■ 6r 2 . 

' r 2 ->r 2 +i5r 2 — Wj — >r 2 

By Stokes' theorem, we can carry out the loop integral over a triangular area encircled by the loop and we have 



(50) 



r 2 — >r2+or 2 — >ri— >r 2 



A(s')-ds' = f 

J a 

~ B(r) 



area encircled by loop 

[Sr 2 x s] 1 



V x A • [Sr 2 x ds'] 



-6v 2 ■ [s x B(r)] 



(51) 



Here a Taylor expansion of the B(r + s') field B(r + s') in powers of s' in the loop integral will lead to terms in power 
of h [HI, [44|. Taking the B(r) field to be located at r in Eq. (|5"Tj) represents the lowest-order h — > approximation in 
the expansion of B(r + s'). As 5 r2 cf) = V T2 4> ■ ^ r 2, we obtain 



V r2 e 



-ie„ f* 1 A(s')-c£s' 



-ie„ /£ A(»')., 



i { e v A(r 2 ) - -[s x e„B(r)] 



and 



Therefore, we obtain 



e, y A(r 2 ) e 



— te„ J r rl A(s')-<is 



= — ie„ f* 1 A(s')-ds' 



-[s x e„B(r)] 



(52) 



(53) 



(-V r2 - e„A(r 2 )) 2 - {( A Vri - e„A(r 2 )) 2 

i —i 



,- ie « G A(s')-ds' i 



je-^/r'a A(s') ds'| 



^ 1 ~, \ 

-Vr 2 - 77 s x e„B r 

i 2 



— v ri + 77 s x e„B(r 

— 2 2 



7,(r,p) 



(54) 
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After some manipulation for the remaining terms, we get 



/ ^ :Fie ~ lP S { mMr ' P)+1 E' Vr/,/(r ' P) + l{ ^ rV + E X ei/B{r)) ■ Vp/,y(r ' P) } = °' (55) 



which is satisfied if 



dtf»(r, P) + | ' V r /„(r, p) + (-V r y + | x e„B(v)) ■ V p /„(r, p) = 0. (56) 



This is just the equation of motion for the distribution function in kinetic theory for a collisionless fluid in weak fields 

Equation (|56p is the equation of motion for the Wigner function of a single-particle state v in external scalar and 
gauge fields. For a many-body system with different states a and particle types v, the Wigner function is 

/(r, P) = n <^f<^( r > P)' ( 57 ) 

au 

The equation of motion for the total distribution function /(r, p) in the collisionless limit will be in the same form 
as Eq. ([5S]) with /„(r, p) in Eq. (|55|) replaced by /(r,p). The presence of other particles in different single-particle 
states allows one use the mean-fields as the external fields and to introduce the collision term by considering residual 
interactions between particles in different av states. 

The exercise in this Section indicates that although the quantum probability density approach is particularly 
appropriate for cases of strongly interacting quantum systems, it also has a well-defined classical weak-coupling limit 
that coincides with the kinetic theory, from which hydrodynamics equations can also be formulated. 



VI. SUMMARY AND DISCUSSIONS 



We have generalized the formulation of Feshbach and Villars to write the Klein-Gordon equation as a set of two 
coupled time-dependent Schrodingcr equations, for the particle and antiparticlcs components of the wave function. 
We have improved upon the formulation of Feshbach and Villars in this re-examination. In our new set of coupled 
time-dependent Schrodinger equations (O the particle and the antiparticle components are better separated, and the 
kinetic energy operator in the equation for the antiparticle component is in the proper form of (p — eA.) 2 /2(E — e v A°), 
with a positive sign. The relativistic properties of the momentum variables in the coupled time-dependent Schrodinger 
equations are more apparent and their connection to relativistic hydrodynamics can be easily established. 

We introduce amplitude and phase functions to cast the time-dependent Schrodingcr equations into hydrodynamical 
form. We find that the equation of motion of the probability fluid of Klein-Gordon particles or antiparticles can be 
written in the form of relativistic hydrodynamics, with an additional quantum stress tensor. The other components 
of the hydrodynamical equation have their classical counterparts. 

For simplicity, we have suppressed the pair-production degree of freedom in the present investigation. The pair 
production is however a quantum phenomenon that can be studied by using the Klein-Gordon equation, as was carried 
out in Ref. j29j to examine the Schwinger mechanism. The presence of the pair-production mechanism is a new feature 
in the hydrodynamical evolution. Since energy and momentum is diverted into pair production, the pair-production 
corresponds to a dissipative process, and will contribute to the viscosity of the fluid. Future investigations to include 
this pair-production in relativistic hydrodynamics will be of great interest. 

As both the Schrodinger equation and the Klein-Gordon equation can be cast into a hydrodynamical form, one 
may inquire whether the Dirac equation can also be written in hydrodynamical form. It is well known that the Dirac 
equation can be reduced to a Klein-Gordon equation, with additional terms involving the spin and particle-antiparticle 
degrees of freedom. For a Dirac particle in an external field we have 

{-f(id„ - eA v ) - (to + S)}^ = 0. (58) 

Upon multiplying this on the left by ^(idy — eA v ) + (to + 5), we obtain 

{(id,, - eA„) 2 - (to + S) 2 - ia ■ eE + a ■ eB(r) - i[-y v d v S]}ip = 0, (59) 

which is the Klein-Gordon equation with additional interactions. Thus, the Dirac equation can be likewise cast into 
a hydrodynamical form, following the procedures outlined in the present manuscript. 

We have thus far discussed particles and antiparticles with opposite charges interacting in a gauge field. For the 
case of neutral particles, there is no interaction with the gauge field characterized by the charges e±. However, 
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rclativistic doubling of states occurs and there are particles and antiparticlcs. How these neutral particles should 
be treated will depend on their interaction with the external field. If the interactions of the neutral particle and 
antiparticle with the external fields arc identical, it will not be necessary to distinguish between a particle and an 
antiparticle. It then becomes possible to construct a simplified theory in which only one degree of freedom enters 
(say, the particle's), with a neutral antiparticle taken to be identical to its corresponding neutral particle. On the 
other hand, if their interactions with the external fields are different, the two degrees of freedom are distinct. One can 
then introduce additional quantum numbers to distinguish particles and antiparticles, and the present investigation 
containing different charges (or quantum numbers) for the new types of interaction will apply. 

In the dynamics of a single-particle state in an external field, the external field can come internally from the single- 
particle state probability density as in the Gross-Pitaevskii equation [3ll - l34| or from the self-consistent mean-field in a 
many-body system [l4|, EU • A many-body system consists of a collection of various particles in their individual single- 
particle states. For example, we can examine a many-body system in the time-dependent Hartree approximation in 
which particle and antiparticles move in a time-dependent self-consistent mean-field generated by all other particles. 
The motion of each particle can be individually followed, as in the time-dependent mean-field approximation in a 
nuclear system [35| . The occupation numbers of these single- particle states will remain unchanged, if there are no 
additional residual interaction between the single particles. When particle residual interactions are allowed as in the 
extended time-dependent mean- field approximation (36j . the occupation numbers will change and will approach an 
equilibrium distribution. As the fluid density is cumulative in nature the dynamics of the many-body system will 
rely on the cumulative effects from particles in individual states. Thus, there will be a total quantum stress tensor 
that is an important part of the equation of state of the many-body system. There will also be contributions from 
the mean fields to the equation of states of the relativistic fluid as in the nuclear fluid [ToT ] . The present investigation 
of Klein-Gordon single-particle states should be useful in the study of the dynamics of a many-body system using 
rclativistic hydrodynamics. 

The foundation of hydrodynamics is usually presented within the framework of kinetic theory, we can compare the 
present quantum probability density approach to the usual kinetic theory approach by determining the equation of 
motion for the Wigncr function using the time-dependent Schrodingcr equation of the Klein-Gordon equation. Our 
comparison indicates that the equation of motion for the distribution function in an external field in the kinetic theory 
can be obtained from the quantum probability density approach in the limit of h — > and weak fields. Thus, the 
quantum probability approach extends the range of hydrodynamics applications to situations where quantum effects 
and/or strong interactions are important. 
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